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Abstract. We investigate the reductions of dispersionless Harry 
Dym hierarchy to systems of finitely many partial differential equa- 
tions. These equations must satisfy the compatibility condition 
and they are diagonalizable and semi-Hamiltonian. By imposing 
a further constraint, the compatibility is reduced to a system of 
algebraic equations, whose solutions are described. 
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1. Introduction 

Dispersionless integrable equations arise in various contexts and have 
attracted people to investigate it from different point of view, such 
as topological field theory (WD VV equation) [331 13 113 ED!, ma- 
trix models [HIEI, conformal maps and interface dynamics fr^ W)\ I41j. 
Einstein- Weyl space (TBI Ej. The hydrodynamic reductions are the 
most developed method to find the exact solutions of the dispersion- 
less integrable equations 0|EJ[I2injl221l2ElllIlllI]- From the hydro- 
dynamic reductions, one can construct the Riemann's invariants and 
the corresponding characteristic speeds satisfy the semi-Hamiltonian 
property or Tsarev's condition [38J. Hence the generalized hodograph 
method can be used to find the exact solutions. Also, the solutions 

l 
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of dispersionless integrable equations can be found by nonlinear Bel- 
trami equation [213 121]; slightly different from the generalized hodo- 
graph method. 

Let's recall dispersionless non-standard Lax hierarchy [33J EH|- Sup- 
pose that A is an algebra of Laurent series of the form 

N 

A = {A\A = aiP% 

t=— oo 

with coefficients dj depending on an infinite set of variables t\ = x, t 2 , £3, • • 

We can define a Lie-Bracket associated with A as follows: 

fA D1 dAdB dAdB . x 

{A, B} = —— - A,BeX 
ox op op ox 

which can be regarded as the Poisson bracket defined in the 2-dimensional 
phase space (x,p). The algebra A can be decomposed into the Lie sub- 
algebras as 

A = A> fc ©A <fc , (fc = 0,l,2.) 

where 



A> fc = {Ae A\A = J2 a iP i } 

i>k 

\ <k = {Ae\\A = J2a iP i }. 



i<k 

Based on this, the Lax formulation of the dispersionless integrable hi- 
erarchy can be formally defined as 

• For k — 0, it's called dispersionless KP hierarchy(dKP) (HH] 

• For k — 1, it's called dispersionless modified KP(dmKP) hier- 
archy [TH El 

• For k = 2, it's called dispersionless Harry Dym (dDym) hierar- 
chy [TTJOEl]- It's the purpose of this paper. 

We define the dispersionless Harry Dym (dDym)hydrodynamic systems 
as follows. The Lax operator of dDym has the form 

X(p) = A_ lP + A + A^ 1 + A 2 p~ 2 + A 3 p~ 3 + ■■■ . 

Then the dDym hydrodynamic system is fH\ El] 

(9A 

(1) — = {X,Q n (p)}, 

where 

n n (p) = (\(pD>2. 
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Here ()>2 denote the projection of Laurent series onto a linear combi- 
nation of X(p) n with n > 2. From the zero curvature equation (t 3 = t 
and t 2 = y) 

~m dy~ = ^(P)M P )}, 

where 

Q 2 (p) = A\p 2 and fl 3 (p) = A 3 _ lP 3 + SA^Aop 2 , 
one can get the dispersionless Harry Dym Equation jTTj 



M_! 3 1 



|2 a-1 / A-ly 



A\ 



dt 4 A_! 
Now, considering the y-flow, we have 
(2) X y = {A, A\p 2 } = 2pA\\ x - (A\) x p 2 \ p , 



or 



A—iy — 2,A_^A 



O.r 



A 0y = 2A± 1 A lx + (At^xAx 



(3) A ny — 2A_ 1 A n +i !X + (n + l)(A 2 _ 1 ) x A n+1 , 

where n = — 1, 0, 1, 2, 3 • • • . Comparing it with the Benney moment 
chain |2*3*1 13ti| IHj. one calls (jHJ) the dDym moment chain. It's the sub- 
ject of this paper. 

The paper is organized as follows. In the next section, one considers 
the reduction problems and obtains finitely many partial differential 
equations. In section 3, we prove the semi-Hamiltonian property when 
using Riemann invariants. In section 4, by imposing a further con- 
straint, one gets some particular reductions. In the final section, we 
discuss some problems to be investigated. 



2. The compatibility conditions 

In this section, we consider the hydrodynamic reduction problems 
following |TJ |ini 122] • For non- hydrodynamic reductions, one refers to 
0. 

We assume that the moments Ai are functions of only iV independent 
variables U{. If the Ai satisfy @, then it's straightforward to show that 
the mapping 

(u-i, u , Mi, u 2 , ■ ■ ■ , Wjv-2) -» (-4-1, A , Ai, A 2 , ■ ■ ■ , A N _ 2 ) 
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is non-degenerate. Hence without loss of generality we set u_i = 
A_i,u = A ,ui = Ai,- ■■ ,un-2 = A N _ 2 - The first iV-moments are 
the independent variables, while the higher moments are functions of 
them, i.e., 

Ak = A k (A_ 1 , A , - • • , A N _ 2 ), k>N-l. 
The equations of motion for A-i, A , • • • , A^-2 become 
dA ■ 

(4) = 2A 2 _ 1 A j+1>x +(j + l)(A 2 _ 1 ) x A j+1 , j<N-3 
r)A 

( 5 ) = 2A\A N ^ X + {N -1){A\) X A N ^ 

oy 

while each higher moment (Ajv-i, • • • ,) must satisfy the overdeter- 
mined system (k > N — 1) using (jlj)© 

~dj = E^;[ 2 ^i^m + (j + i)(^-i)^ + i] 



0A k [2A2 JfdA^ d A^ + {N _ MxA 



N-l\ 



dA N _ 2 ^ dAj dx 

= 2^_ 1 [x;^ ±i ^]+(*+i)(^ 1 )^ +1 . 

j=-i 'j 

Comparing the coefficients of ^f(j = —1, 0, 1, • • • ,N — 2), one has 

, . dA k+1 dAu dAu dA N ^ 

(6) — *±i = *_ -| 5 ?-L o < i < N - 2 

1 J dAj dA^ dA N - 2 dAj ' " J ~ 

(7) ^aZT = ^dA^ + l)A ^ + A ^ dA N .2 d~A~! 

- (A; + l)A fc+1 
Now, letting k = N — 1 and defining r = logA_i, one has 

{ ] dAj dAj^ dA N _ 2 dAj ' - 7 " 

(9) = £ -a^-o + + 9An _ 2 dr 



NA 



N 
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The compatibility of (jHJ and (jHJ) gives a system T of Af ^~ 1 - > non-linear 
second order equation for the single known Ajv_i(A_i, A , A 1; • • ■ , v4jv_ 2 ) 
One can show that by induction if T is satisfied then the analogous 
compatibility for A^ik > N) is also derived. Let's investigate the case 
N = 2 in more details. Then we have (A = s) 



ds 
dA 2 

dr 



A lr exp(-r)+A( s 
A X A U + A ls A lr - 2A 2 . 



Cross-differentiating, one gets the quasi-linear second differential equa- 
tion 

(10) Ai„ exp(-r) + A ls A lrs - {A x + A lrs )A lss + A lr exp(-r) 
+ (A ls ) 2 = 0. 

Letting 

A\ = a, A lr = b, Au = c, 

one can express (fTUj) as the degenerate non-homogeneous hydrodynamic 
system 










\ 







— cexp r 


(a + b) exp r 


! 





1 








Define the characteristic speeds as 



u 



-c exp r + a/ c 2 exp(2r) + 4 (a + b) exp - 



-cexpr + A/c 2 exp(2r) — 4(a + 6) expr 



A simple calculation yields, 
(12) 



u \ v OWn 

u / I u I \ v 



v—u 

UV I 5 



there being no a-term! It is of non-homogeneous hydrodynamic systems 
of Tsarev-Gibbons type |2"T) ITH] and it has one obvious hydrodynamic 
type conserved density (u + v) . Moreover, according to the theory of 
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Poisson commuting Hamiltonians [!%] . one can also find a conserved 
density of first derivatives: 



For a given affinor v\ [u] of a hydrodynamic system 

< 13 > ^ = <^% 

one can define the Nijenhuis tensor N^Au) 

A7fc .^i fc^j 

« 4 du s 3 8u s s d Uj s dm 
Then we can define the corresponding Haantjes tensor H l - k {u) [21] 

H} k = (nix - K v >" ~ ( N >1 - K 3 <H- 

The system f|13|) is diagonalized if and only if the Haantjes tensor 
Hj k (u) vanishes identically and all the eigenvalues of the affinor vj(u) 
are real and distinct. That is, there exist iV-functions A n (Riemann in- 
variants), depending on the variables Ui, in which the equation (fT3|) is 
diagonalized 

(14) ?± = vJ K 



oy OX 

where V n are the eigenvalues of the matrix v\ (u) , called characteristic 
speed. For the hydrodynamic system (j3J) and (0) of the reduced dDym, 
the Haantjes tensor H l - k {u) vanishes identically whenever the system 
©and (jnj) is satisfied. Therefore, any consistent reduction of dDym 
is diagonalizable and then we can use Riemann invariants to discuss 
problem further. 

Finally, one remarks that as in the case of dKP hierarchy [2H], a 
similar argument shows that for the reduced dDym hierarchy (JTJ we 
also have the Kodama-Gibbons formulation: the Riemann invariants 
are 

d\ 

(15) A n = X(q n ), where — (q n ) = 0, n = 1, 2, - • • , N 

op 

and then the hierarchy (JTJ can be expressed as 

i)K n m {v n ) 9K ' 



where V n = (q u q 2 ,--- , q N ) and Q m (V n ) - 



a. " "III \'Tll r\ 

Ot m ox 

q 

dp \p=V n ' 
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3. semi-Hamiltonian property 

In this section, one will prove the semi-Hamiltonian property (Tsarev's 
condition) [35] of the reduced dDym hierarchy using Riemann in- 
variants. Suppose that each moments A n can be expressed as the Rie- 
mann's invariants Aj, which satisfy the equation(t 2 = y) 

(16) — = V — 

dy % dx ' 

where 

Vi = 2A\q t , i — 1, 2, 3 • • • ,N. 
Then the moment equations © can be written as 

(17) ViA nX = 2A\A n+lx + (n + l)A n+1 {A\) Xi , 

On the other hand, from ((TJ), we also have, equivalent to (fTTj). 

dX _ 9X,^2 \ p 2 2 ^ A 2 p 
dy dp ~ l x ^ dx 

and then, using ()16|) . one obtains, after reshuffling terms, 

9X 2 dX, , 2 , 1 

(18) Wr P lTp iAUh -V^2A^ P 

Cross-differentiating, one has (0 = A 2 _^) 

^(Vi + 2<Pp)(V 3 + 2<Pp) 2 - <p Xt <p X] 2<Pp 2 (y t + 2#) - {Vj + 2<pp) 2 <p Xi 
^<P X] {V j + 2ct>p) 2 



m 

dX j 
Qy. 

= ct> XiX] {V 3 + 20p)(^ + 2<p P ) 2 - Ai Aj 20p 2 (^. + 2<Pp) - (Vi + 20p) 2 0A^ 

-2p<p Xi <p X] (y i + 2<p P ) 2 . 

Letting p = we obtain 

On the other hand, comparing the coefficients of p-power and using 
(TTTIj) . we get the only equation 

(20) ^ = —^JV^W 2+1] > 

The higher moments A n , with n > can be solved recursively using 
(JTTj). These equations (fE?j) (}2Tjj) are compatible and their solutions are 
parameterized by 2N functions of a single variable. 
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A direct calculation, using MAPLE, confirms that the reduced equa- 
tion is semi-Hamiltonian, that is, 




2 V V K-K: 



9 I d\ k 



OX, \ Vj-V^ 

for i,j, k all distinct. Then the reduced equations (fTBj) are thus inte- 
grable by the generalized hodograph transformation 



4. Algebraic Equations and Special Reductions 

To investigate the reduction problems, we introduce A = 2\nA^i to 
put (fTT?|) and (}2T)j) in a more compact form(z ^ j) 

dv > l A r W > +1/1 
A XiXi = A Xi A Xi [,Xj 



or, noting that Vi = 2A 2 _ 1 q i , 
(22) A XiX] = A Xi A 



QiQj 



3 (Qj ~ 9i) 2 ' 

Then as in [21] for the dKP case, we impose two further restrictions on 
the reduced system (j2"Tj) and (|2^|) . First, from the form (|2*Tj). we require 
that the reduced system is translation-invariant in the sense that 

A; — > Xi + c => qt —> g i5 A -> A 

or,equivalently, 

(23) = 

(24) cL4 = 0, 

where 5 = YliLi ^fr- Secondly, we require the homogeneity of the 
functions A and qj in the variables Aj. A should be of weight 0, and 
the of weight —1. Hence 

(25) Rq t = —q { 

K 

dA dA 

(26 » V = -ox- 
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where k is a positive integer and R = Y^i=iXi-^-- Plugging (J2~2~j) into 
(f!£Ij) ()26)) and eliminating the second derivative, we obtain 

dA OA A d 2 A d 2 A 



3=1,3^1 J 1 

^ j dX i X j + { ^ dXiXj ' 
3=1,3+% J j'=1j¥* 

N 

Hence, either A\ ( = or 

TV 

Similarly, plugging (}2*T]) into (J2~3~)) and (J2SJ), we get 

TV Q TV Q 



9A< 



3=1,3+%' 3 



(28) 



1 m \ ^ ^ r ft^i 

2?, (VA,) a- | «"'' 



The systems (|2*7jl and (J28|) form a 2iV algebraic equations for 2iV un- 
knowns, the qi and J-^. The solutions of the system are essentially 
unique. With two Riemann invariants and k — 1, a simple calculation 
can yield by carefully choosing the integration constants 

(29) „ = -,(, = _« ^ 1 = ^i_M. 

Ai — A2 lo 

Hence from (fT6|) we have 
Ai 



(30) 



32 

A 2 /., V 






(x,\ 


)l 


\ A2 / 
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This equation (|HUJ) will correspond to the reduction of Lax operator 
constructed by Riemann-Hilbert method fH\ IT2] 

(31) A(p) = A_ lP + A + -. 

V 

If we consider more general Lax reductions of the form for positive 
integers m and k [TTJ [T3] 

(32) \(p) m = A™ 1 p m + w m _ lP m - 1 + w m _ 2 p m - 2 + --- 

+ W- K+1 p^ +1 + p~ K , m + K = N, 

then from it's not difficult to see that (|32|) satisfies the conditions 
(J25J) ({21} (j25j) (j2Sj) ■ Hence its corresponding characteristic speeds and 
A_\ will be one of the solution of the 2N algebraic equation (J27|) and 
(|2*5|) . It would be interesting to know whether any other solutions exist. 

Next, one generalizes the operator R in (|25|) or (|26|) to the following 
forms 

N 

R = ^2g k (\)d Xk h k $) 

k=l 

(33) RVi = —Vi 

K 

and 

N 



(U) R M - 

(34) R dx, ~ ~5A? 

where A = (Ai, A 2 , • • • , Ajv) and (g k , h k ), (g k , h k ) are arbitrary functions 
satisfying the conditions compatible with 



(35) Sg k = Sh k = 5g k = 5h k = 0. 

Then a similar calculation can get the following 2N algebraic equations 
generalizing fl2Z|) flU 

N nn 1 N Bh 

(36) £ (gA - &fc)^[-ML- + ft] = _ + «J>-!) 

3=1,5^1 IJ .7 = 1 J 

(37) £ Q?A- - <M^t-^ = -(i + E^)- 

For N = 2 and k = 1, simple calculations can show that a suitable 
choice of (gi, hi) , (g 2 , h 2 ) , (gi, hi), (g 2 , h 2 ), not unique, can also obtain 
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the solution Iffijl . 

Finally, one notices that if we let 

hj = h 3 = l 

for j = 1, 2, • - • , N, then we get the weaker condition than (|33j) 

Sgk = $9k- 

Hence the equations (|3^j)(j3Tj) reduce to 

£ a = 

Hi ~ Hi K 

2V 

V (pj - gi)A\ gtgj = -1 

If = gj = Xj, then we can obtain ()27j) and (|2*Hj) . 

5. CONCLUDING REMARKS 

We prove the semi-Hamiltonian property of reductions for the d- 
Dym and find some particular solutions invariant under translation and 
homogenity. In spite of the results obtained , there are some interesting 
issues deserving investigations. 

• The integrability and solution structure of the equation (fT^j) (or 
(fTUm is unclear dl EH • It is not difficult to see that (JT2J) 
can be extended to : 

dui dui uiu 2 • • • it 7v , / N 

^- = n— + - — -, where n = ( ) u k ) - Uf. 

dr ds Ll k ^i(uk - Ui) 

One hopes issue these problems elsewhere. 

• In jsJI], the algebraic reductions for dKP are found and in 
[23 ESI, the waterbag reduction(non-algebraic) for dKP is also 
found. If we define 

/oo 
q n f{q,x,y)dq, n = -1, 0, 1, 2, • • • 
-oo 

then we obtain 

(38) X=p 2 (P ^-dq), f = f(q,x,y), 

J-oo P-Q 

where P J denotes the Cauchy principal value of the integral. 
Also, from (|HS|). the "distribution" function f(x,y,q) must sat- 
isfy the Vlasov-like equation [21 E3 El 

(39) f y = 2A\qf x - {A\) x q 2 f q = {/, A\q 2 } x , q . 



12 



JEN-HSU CHANG 



Comparing and (j2J), we can assume / = F(X) for any 
function F. Hence the Lax operator A will satisfy the non- 
linear singular integral equation 



This equation will help us find the (no n-) algebraic reductions 
[221 EH], study the initial value problem for d-Dym as in the 
case of dKP (121 US] and it needs further investigations. 
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